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1. Introduction 

20 years after the investigation of transport in stratified random media by Matheron 
and de Marsily [|I], their model is still intriguing. Not only because it is exactly solvable, 
but also as it serves as a prototype of models exhibiting super-diffusive behaviour 0. 
As such it found applications in various fields like hydromechanics ||, f§] and polymer 
physics i, i. 

The original model can be described as follows: Consider a two-dimensional random 
medium consisting of distinct parallel layers (or strata) which extend in xi-direction 
and have varying permeability. The resulting velocity field u for a fluid flow in this 
medium depends only on the component £2 transverse to the layers. In a random 
ensemble of such stratified media, it is a random field with given statistical properties. 
The properties of such two-dimensional stratified models are well known in the 
asymptotic limit of large times, t — *■ 00, and for the case of vanishing local diffusion, 
see e.g. 0, ||, ||. The large scale transport behaviour is roughly characterised by 
the mean square displacement of the tracer particles, appropriately averaged over the 
ensemble of the random flow fields. The time derivative of this quantity encodes the 
large scale diffusion or, as it is called, 'dispersion' coefficient D(t) of the tracer cloud 
as a function of time. The transport is called diffusive if this coefficient becomes a 
finite constant in the asymptotic long-time limit t — > 00, it is called super-diffusive if it 
asymptotically grows with time as D(t) ~ t a with a > 0. For a flow field parallel to the 
strata the model exhibits (for t — ► 00) such a super-diffusive behaviour with a = 1/2. In 
the more general case, where one also has a transversal flow component, the dispersion 
coefficient eventually tends to a constant value, i.e. the transport shows normal diffusive 
behaviour. 

However, what is still lacking, to the best of our knowledge, is a detailed analysis 
of the pre-asymptotic behaviour of the model and, even more interestingly, of the 
corresponding sample to sample fluctuations of the dispersion coefficient. The latter 
is crucial for the reliability of the model. Small fluctuations indicate that the transport 
parameters which are constructed as averages over the random ensemble are indeed also 
representative for the behaviour of individual realizations of the ensemble. 
In the present article, we try to close this gap. Assuming a stratified Gaussian random 
flow with a Gaussian correlation function, we derive exact expressions for the dispersion 
coefficient and its mean square fluctuations as a function of time. The latter quantity 
encodes the sample to sample fluctuations of the dispersion coefficient. We show that it 
remains finite even in the asymptotic limit of infinite times t — > 00 as long as the initial 
tracer distribution is of finite extent which implies that in the given two-dimensional 
model the dispersion coefficient, surprisingly, is not self-averaging. 
The paper is organised as follows: In Section ^| we give a brief description of the 
two-dimensional Matheron-de Marsily model and define the corresponding transport 
parameters. The next two sections deal with the exactly stratified model (i.e. a model 
where the flow field is aligned parallel to the strata), and with the more general case of 



Stratified disordered media 



3 



a non- vanishing transverse flow component. In Section |5] we analyse the self- averaging 
properties of the given transport parameters. 

2. The model 

The transport of a passive tracer particle by a two-dimensional flow field u(x) is modelled 
by a Langevin equation dx(t)/d£ = u(x) + where x(i) denotes the position- vector 
of the particle, and £(t) is a Gaussian white noise which generates the local diffusion 
process. It has zero mean (&(£)) = (with % = 1,2), and a correlation function 
(£i = 2Dij6(t — t'), where (...) stands for the average over the white noise 
ensemble. We assume a diagonal local diffusion tensor, = 8{j Di which, however, 
may have different entries for the longitudinal and the transversal diffusion constants, 
D l + D 2 . 

In a two-dimensional stratified medium the flow field u depends on x 2 only. In the 
model investigated here, one has u(r) = (u + u'(x 2 ),v Q ), where u and v denote 
the components of a constant drift velocity in 1-and 2-direction, whereas u'{x 2 ) is the 
random flow contribution generated by the layer to layer permeability variations (see 

e.g. I). 

The Langevin equation therefore reduces to 

j t x 1 {t) = u +v!{x 2 {t))+tx{t) (2.1) 

j t x 2 (t) = v + &(t) . (2.2) 

The flow contribution u'(x 2 (t)) is a random function chosen from an appropriate 
stationary ensemble which defines its stochastic properties. In the model investigated 
here, we assume that it has zero mean and a Gaussian correlation function, 



u'(x 2 ) = (2.3) 



(x 2 — X 



I \2 S 



u'(x 2 )u'(x' 2 ) = ai exp \- ^ 2p l > j (2.4) 

The over-bar denotes the average with respect to disorder ensemble, the parameter 
cr quantifies the disorder strength, and / denotes the correlation length of the field. 
As is shown below, these two relations uniquely determine the temporal behaviour of 
the dispersion coefficient, higher correlation functions do not enter. This, however, is 
no longer true if one deals with the sample to sample fluctuations of this dispersion 
coefficient. In this case one needs also an explicit expression for the four-point function 
of the random field. For simplicity, and without severe restriction of generality, let us 
assume a Gaussian ensemble, so all higher correlation functions can be decomposed into 
products of the given two-point function. 

Equation ( |2.1fl and ( |2.2j ) define a system of differential equations which is readily solved 
by 

*i(*)|xo = x 01 +u t+ f &t'u\x 2 {t')) + [ dt' Ci(t') (2.5) 

Jo Jo 
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X 2 (t) |xo =Xo2+Vot+ I dt' £ 2 (t') • (2.6) 

Jo 

The single tracer particle subject to a given noise history is characterized by its 
trajectory x(t) | Xo , where the subscript x indicates the fact that the tracer coordinates at 
a given time t still depend on the initial injection position x(t = 0)| Xo = x = (x i, ^02)- 
In the general case where one has a spatially extended injection of material at time t = 
the injection points x are distributed according to a distribution p(x ) which without 
restriction is chosen to be normalised, 

J d 2 x p(x ) = 1. (2.7) 

The shape of the resulting concentration distribution is characterised by its moments 

(Xi(t)») = Jd 2 x oP (x ) (xiWlS,) (2.8) 

(with n = 1,2,...) which in particular are used to define the centre-of-mass velocity 
and the dispersion coefficients of the resulting concentration plume. The centre-of-mass 
velocity (for a given realization of the flow field) reads 

Ui(t) = j t (xi(t)) . (2.9) 
The dispersion coefficients which give the spreading rate of the plume are defined by 

Dij(t) = \^ t { fa® X M) - fo(*)> (*;(*)> } (2-10) 

In a stochastic approach, the given single medium is viewed as one particular realization 
of a spatial stochastic process. Large scale transport coefficients are derived from 
appropriately constructed averages over the ensemble of all possible medium realizations. 
These averages represent statistical properties of the (artificial) ensemble as a whole. 
At first glance, they therefore seem to be of limited predictive value with respect to 
the properties of a single given aquifer. For appropriately chosen quantities, however, 
in many cases the fluctuations from realization to realization become small as soon 
as the cloud has sampled a sufficiently large representative part of the given medium. 
As soon as the transport properties found in different realizations of the medium then 
fluctuate only weakly around the corresponding ensemble averages, these averages again 
also predict properties characteristic of a single typical realization of the medium. 
In what follows, we investigate such ensemble averages of the centre-of-mass velocity 
and the dispersion coefficients. We call the averaged quantities corresponding to ( |2.9| ) 
and ( p,10| ) "effective" quantities. The resulting effective centre-of-mass velocity and 
effective dispersion coefficient are defined as 



and 



uf (t) = ui(t) = - (x t (t)) (2.11) 

Df{t) = D~iT) = ^{Mf - (Xi(t)) 2 } . (2.12) 

The non-diagonal dispersion coefficients vanish due to symmetry reasons. Note that 
two types of averages are involved, the average over the white noise which generates 
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the local diffusion process, indicated by the brackets, and the average over the disorder 
ensemble, indicated by the over-bar. The order in which these averages are performed is 
crucial, as the following slightly different definition of an averaged dispersion coefficient 
demonstrates. We define what we call "ensemble dispersion coefficients" as 

DT(t) = \^ t {T^-J^W) 2 } • (2.13) 

This quantity, which also is well known and frequently used in the literature, represents 
the dispersion characteristics of the whole ensemble of realizations of the flow field. It 
takes into account an artificial dispersion effect caused by fluctuations of the centre-of- 
mass positions of the tracer distributions in different realizations of the inhomogeneous 
medium. This effect is suppressed in the effective dispersion coefficient as defined in 
( |2.12| ) because there the centre-of-mass positions are superimposed before the ensemble 
average is performed. For an extensive discussion of these problems and the related 
literature, see [|TU]| and [II |. In general, the experimentally observable dispersion which 
is a property related to a single given realization is represented by the effective quantity 
D eS (t). Usually, at finite times the ensemble quantity overestimates the experimentally 
observable dispersion considerably, see again the detailed discussion given in fTOl , fTT|[ . 
In the asymptotic long-time limit t — > oo, however, one usually expects that the two 
types of dispersion coefficients become equal. 

For the model investigated here, the full temporal behaviour of the given transport 
coefficients can be evaluated explicitly. The corresponding calculations are elementary 
but somewhat tedious. They are summarized in |Appendix A| . Using the explicit solution 
( |2.6| ) and performing the appropriate averages, one finds the more or less trivial results 

uf(t) = u , uf{t) = v and D™(t) = D*{t) = D 2 (2.14) 

valid for all times. A more interesting behaviour follows for the longitudinal dispersion 

coefficients. Here one derives the relations 

1 



DlT(t) = D, + -d t I I u{k)u{k>) x 



k k 

rt 



X 



/ dt' / ^t" e ~( k2D 2i'+k' 2 D2t"+2kk'D 2 rain(t',t")) 
JO JO 

DfAt) =D™(t)-U t Uu(k) Jd2x oP (x )e- ik y° J*dt>e~ k2D ^ 2 (2.16) 

where we use J k . . . = (27r)~ 2 / d 2 k ... as a convenient short-hand notation. 
One can already observe three simple features: First of all, the effective and the ensemble 
dispersion coefficient in the direction transversal to the flow field are both equal to the 
local dispersion coefficient D 2 . The heterogenities of the model do not change this 
property because they do not depend on the fluctuating part of the flow field. The 
particles perform a simple random walk in 22-direction. Secondly, the effective dispersion 
coefficient in xi-direction does not depend on the longitudinal extension of the initial 
distribution. This stems from the infinite correlation length of the flow field in that 
direction. The third observation is that the ensemble dispersion coefficient does not 
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show any dependence on the size of the initial distribution. 

What remains to be evaluated is the explicit behaviour of the longitudinal dispersion 
coefficient. In the following we will discuss the cases vq = 0, where the fluid flow is 
aligned with the strata of the medium (the so called 'exactly stratified model'), and 
«o^0 separately as they give rise to a different asymptotic behaviour. We will further 
distinguish between a point-like initial distribution at x 10 = x 2 o = (which implies 
p(x ) = 5 (x ) in (|2.8| )) and an spatially extended distribution 



P(Xo) = 72¥L eXP l"f§j 5(X0l) (2 ' 17) 



where the injection points are aligned along a line in x 2 direction. As it turns out, the 
line shape is no restriction of generality since a finite extent in xi-direction does not 
alter the given results. 

3. The exactly stratified medium 

For v o = and a point-like initial distribution the expressions for the effective dispersion 
coefficients given by equations (|2.15|) and ( |2.16| ) can be evaluated easily. One derives 



Dl?(t) = D 1 + c- 2 T D x 1 + --1 (3.1) 




Dg{t) = D™(t) - a\ r D ( + ^ - ^1 + ^ J (3.2) 
where the time-scale td = I 2 / D 2 is a measure for the diffusive transport over one 



transversal correlation length. The result is plotted in Figure [Bl]. There are two 
characteristic regimes. For t the dispersion coefficients asymptotically are given 

as 



DlT(t) = D 1 + a z Q t + O {(t/r D Y) (3.3) 

D«{t) =D 1 + 0((t/r D ) 2 ) (3.4) 

In this regime there exist great fluctuations between the centre of mass coordinates of the 
source in different realizations of the flow field, the transport is dominated by advection 
which explains the linear time dependence of the ensemble dispersion coefficient. In 
each realization the particles only experience local diffusion which can be seen in the 
effective value which is approximately constant. This behaviour characterises the Taylor 
regime as already described by [[§. 

In the opposite long time limit, t 3> t d , the dispersion coefficients take the following 
asymptotic form 

Dtf (t) = Di + a 2 Q t 1 / 2 + O ((t/r D y 1/2 ) (3.5) 

B$(t) =D 1 + a 2 [2 - y/2] t 1 ' 2 + O ((t/r D )- 1/2 ) (3.6) 

They show the well known anomalous diffusion t 1//2 -behaviour with different pre-factors 
for the two types of dispersion coefficients, a result which is already described in 0, 
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and ||. These authors derive this result by considering the asymptotic behaviour of a 
particle exhibiting a random walk and experiencing a different velocity after each time 
step, i.e. in each layer. The factor of (2 - y/2\ between D<$ s (t) and is due to 

large fluctuations in the centre of mass coordinates of the source distribution in different 
realizations. 

For an extended initial distribution (in particular for the line source given by equation 
( |2. 17] ) the expression for the ensemble coefficient D™ s (t) as given in (|3.1| ) remains 
unchanged. The effective dispersion coefficient, however, now reads 

D$(t)\v nB = D«"(t)- 





2t 



T~D+T L V T D -\-T L 

2 



(3.7) 



Another time-scale that becomes relevant, 17, = L / D2, is the typical time for diffusive 
transport over the width L of the initial concentration distribution. For times much 
smaller than 77 the effect of a finite initial extent of the injection region is equivalent 
to averaging over various realizations. The dispersion coefficient tends to the ensemble 



value given by given by equation (3.1). For times much longer than 77 the system 
behaviour no longer depends on the initial width of the injection region. The effective 
dispersion coefficient, therefore, tends to the value derived for a point source, equation 
Q3.2|) . The resulting temporal behaviour is plotted in Figure [B2| . For an infinite L the 
value of I?if(t)|]ine fends to D™ s (t) for all times. 



4. The general case 



Unlike the exactly stratified model, the long-time behaviour for the case v 7^ shows 
normal diffusion for asymptotically large times, t — > 00. Evaluating equations ( 2.15Q 
and ( 2.16 ) now lead to the more complicated expressions 



Dtm=D 1 + ] l^c- 2 T v erfL 



erf ^ 




l + 2t/r D ) 





+ 



+ exp 



D 



erf 



t2 l 2T " u it+r 



\ 1 + 2t/r D 



D^{t) = D-{t)-^alr v 



+ exp ^(l + 4t/r D ) 



erf 



\ 1 + 2t/r D 



erf 



+ 



T: 



\2 



D_ 

r, 2 



erf 



t 2 K 
\ 1 + 2t/r D 



(W* + 3) - 



— erf 



T 



\2 



r, 2 



1 + At/r D ) 



(4.1) 



(4.2) 



where erf (x) is the Gaussian error function as defined in The time-scale r v = 1/vq is 
set by the time necessary to transport a tracer particle advectively over one correlation 
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length of the flow field. 

For td <C t„ three regimes can be distinguished which is shown in Figure [B~3| . For small 
times the advection has only a negligible influence on the dispersion coefficient. We 
can therefore identify the two superdiffusive regimes from the exactly stratified case, 
i.e. the Taylor-regime given by eqns. ( |3.3| ) and ( |3.4j ) for t <C td "C t v (neglecting terms 
of the order 0(t d /t v ) and O ((t/r^j^resp.) and the superdiffusive t 1//2 -behaviour 
given by eqns. ( |3~5|) and (|3.6|) for t d <^ t <^ t v (neglecting terms of the order O (r D /t) 
and O ((£td/ t ^) 2 ) , respectively). In the asymptotic limit t — > oo the system exhibits 
normal diffusion. The values of the dispersion coefficients are 



D e ™(t) = A + V2ir a 2 r v + (r D /r v ) + O {t v /VWd) (4.3) 

Dfi(t) =D 1 + ^V2^a 2 r v + O (t d /t v ) + O (r v /y/fo) (4.4) 

which shows that, somewhat unexpectedly, the two quantities do not become equal 
even in the normal- diffusive case. This model seems to be a counter-example to 
Metzger's conjecture [I3| that for normally diffusive models the effective and the 
ensemble dispersion coefficients should always become equal in the asymptotic limit. 
One remarkable feature of the cross-over between the superdiffusive and the normal 
diffusive regime is that it does not occur at r v , but at a time scale r^/r D = D 2 /vq which 
is independent of the correlation length of the system. 

For the case t v <C td the behaviour is much simpler. Advection dominates over diffusion, 
i.e. the diffusive time scale td cannot be resolved. Therefore there are only two regimes: 
for t we again find the Taylor regime with a linear growth of the ensemble dispersion 
coefficient, at r v there is a cross-over to a normal diffusive regime given by 

DlT(t) = D 1 + x flalr v + O (^) + O ( ,[^\ (4.5) 



2 

Because the diffusive effect on the mixing is suppressed by advection, the effective 
dispersion coefficient is equal to D\ in leading order: 

D™{t) = Dl + 0( T 4)+0 (J?±\ (4.6) 



5. Self- Averaging 

The question remains as to whether the effective ensemble-averaged parameters 
represent the typical behaviour in a single realization of the statistical ensemble. To 
investigate this problem further, we calculate the root-mean-square sample-to-sample 
fluctuations of the centre-of-mass velocity and the dispersion coefficient which are given 
by 

(5u(t)f =^t^-Mt) 2 (5.1) 
(5D(t)f = "Dffl- ~W) 2 (5.2) 

where u(t) and D(t) stand for the centre-of-mass velocity and the dispersion coefficient 
of the single realizations as defined in (|2.9|) and ( p.lOj ). 
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Again, these quantities can be evaluated explicitly in the given model. As expected, 
we find that the effective velocity becomes self-averaging in the asymptotic long-time 
limit, Su(t) — > for t — > oo. In the opposite short-time case, t — > 0, the fluctuations 
are determined by the fluctuation strength of the disorder <Tq. As an example, consider 
the simple case of an exactly stratified medium, vq = 0, with point-like injection (i.e. 
L — 0). The explicit temporal behaviour of the root-mean-square fluctuations is found 
to be 



(Ht)) 2 (M*)) 2 I 1 f53l 

Mtf V ul \/l + 4t/r D • 1 • ) 

Surprisingly, the situation is different in the case of the dispersion coefficient. It does 
not become self-averaging even in the asymptotic long-time limit. Here, in the limit 

2 

t — > oo the normalized fluctuations 6D/D tend to a constant finite value for the exactly 
stratified case v o = as well as for the case v o 7^ where on has normal diffusion in this 
limit. 

The explicit results are somewhat cumbersome and not very illustrative. The derivation 
and the explicit formulae are, therefore, deferred to Appendix B| . Figure [B4] illustrates 
the qualitative behaviour for the case of an exactly stratified medium (i.e. vq = 0) and 
line injections of different lengths L. As already described in the previous two sections, 
the short time regime is characterized by an effective dispersion coefficient which is 
mainly determined by local diffusion. This quantity does not vary between different 
realizations of the flow field, the fluctuations therefore obviously vanish for t — > 0. For 
t — ► 00 the sample-to-sample fluctuations tend to a constant value, irrespective of the 
initial length of the injection source. The time-scale on which this asymptotic value is 
reached is given by tb + Tl, i-e. it is larger for larger L. Consequently, an infinitely 
long source will show no fluctuations at all as all of the fluctuations of the flow field are 
already experienced right from the start. 

The non-vanishing sample-to-sample fluctuations demonstrate the fact that the effective 
dispersion coefficient as defined by ( |2.12| ) is not a self-averaging quantity, it therefore 
does not represent the typical behaviour expected in a single representation of the flow 
field in the two-dimensional stratified model. 



Appendix A. 

In the following appendix, we sketch the basic steps used to derive the results for the 
temporal behaviour of the transport coefficients discussed in the text. Starting point 
is the solution of the Langevin equation given by ( |2.5|) and Q2.6|) which is inserted into 
equation (|2.8| ) with n = 1 and n = 2. To perform the necessary white- noise averages, it is 
useful to replace the flow field by its Fourier representation, u'(x2(t)) = J k u'(k) e~ lfc:C2 W. 
After inserting ( |2.5| ) and (|2]6|) , the resulting white-noise averages can be calculated using 
the general formula (exp{-i £j $&t'Lj{t') = exp{- Ej,iJd? Lj(t') DjiL^t')} 
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and an appropriate choice of the arbitrary auxiliary function Lj(t). One ends up with 

rt 



jk(x 2+v t') e -k 2 D2t' 



< %i (t) |xo > = x i + u t + / dt' / M / (A;)e 1 

JO Jk 

= 2Dxt + (x i + u t) 2 + 

+ 2(x 01 +u t) f dt' f u>(k)e ik{xo2+vot ' ) e- k2D2t ' 



(A.l) 



<x 2 {t)\ 



xo 



> 



+ 



+ 2 / dt' 



C dt" { { s'(A;)u / (A ; , )e- ifc(xo2+t,0 *' ) - ifc ' (:E02+1 ' '" ) x 

Jk Jk' 

exp (-k 2 D 2 t' - k' 2 D 2 t" - 2kk'D 2 min(t', t" )) (A.2) 



and 



< 2?2(t)|xo > = X 02 + Uo* 
< ^WUo > = 2^2^ + (x 02 + VQtf 



(A.3) 
(A.4) 

These expressions still depend on the initial positions x of the particles. According 
to (2£) one finally performs the remaining average over the initial distribution p(xo). 
Inserted into the definitions for the effective and ensemble dispersion coefficients (|2.12| ) 
and ( |2.13|) , and performing the ensemble average over the disordered flow field, this 
eventually leads to ( 2.14 ), ( p. 15 ), and ( 2.16 ). Using the explicit form for the initial 
distribution of equation (|2.17|) and the expression for the disorder averages ( [2.3|) and 
(|2.4|) then yields the final explicit results given in the text. 



Appendix B. 



In [Appendix B| , we sketch the basic steps used to derive the results for the sample-to- 



sample fluctuations of the transport coefficients given in the text. 

For an exactly stratified medium (i.e. vq = 0) and a point-like injection (i.e. L — 0), 

using the steps described in |Appendix A|, one derives the following results for the square 



of the dispersion coefficient D(t) in a given realization of the medium as defined by 

D{tf = D\ + (B.l) 
+ 2A / / u(h)u(k 2 ) /*dt ie -^ D2 *- fc ^2ti (g-^faUat! _ A + 

Jk! Jk 2 Jo V ' 

+ / / / / u(k 1 )u(k 2 )u(h)u(h) x 

Jk\ Jk2 Jk% Jk4 

x ^-k\D2t~k\D 2 t! (^-2k x k 2 D 2 tx _ ^ e -k%D 2 t-klD 2 t 2 ^ e -2k 3 k 4 D2t2 _ A 

Assuming a Gaussian ensemble for the flow field, the four point function 
u{k\)u{k 2 )u{k^)u{ki) factorizes as usual into products of two-point correlation functions. 
Using this decomposition, and a rescaling of the time integrations of the form T\ : = 
ti/t we get for the root-mean-square sample-to-sample fluctuations of the dispersion 
coefficient defined by ( |5.2|) 

(5D(t)) 2 = \4mt } dn j dr 2 [a ,n,r 2 )~ 1/2 + B g, Tl>Ta )" 1/a 

^ 
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- 2 



A 



2t 



Ti,T 2 -T 2 



-1/2 



+ 



2t 



n, r 2 ) - (ri + r 2 )' 



- 2 

1-1/2 



5 



+ 



Tp 
2t 



B 



• n,r 2 - t 2 



-1/2 



td 
2t 



Ti,r 2 ) - (ti + r 2 ) z 



-1/2* 



(B.2) 



with 



A 
B 



IE 
2t 
td 
2t 



,Tl,T 2 

,n,T 2 



T_D_ 

2t 

TD_ 

2t 



+ 2 



td 
2t 



+ Ti+T 2 



+ 1 + Ti 



2t 



+ 1 + r 2 



(B.3) 



The calculation is easily generalized to the case of an extended injection (i.e. L ^ 0). 
There, the functions A {^-, r l5 t 2 ) and B (j^-, Ti, t 2 ) are replaced by 

,TD + T L 



A 
B 



Tp + T L 

2t 
Tp + T L 

2t 



,Tl,T 2 



■ 2t 
Tp + T L 



+ 2)(^ + r 1 + r 2 ) 
+ l + ri)(^^ + l + r 2 ) 



~ ,Tl,T2 ) :=[ ^2r L + l + Ti)[ ^2r L + 1 + T2) (R4) 

For a medium with a drift component perpendicular to the strata (i.e. v ^ 0) and a 
point-like injection (i.e. L = 0) , we get 



(5D(t)) 2 = J J C(h)C(k 2 ) 



ki ki 



■ e -2k*p 2 t F(k u k2)F*(k u k2) 
+ e -i,o(fei-fc 2 ) ie -(fe?+fc 2 2 )D 2t F(A;i! A; 2 )F*(A; 2 , fci) 
where C(A;) is the partially integrated correlation function 

k 2 l 2 



and 



C{k) :- 
F(h,k 2 ) 



k-^ J J k\ 



u(k)u(k') = v27r<7 /exp 



(B.5) 
(B.6) 



1 



iv k 2 + k 2 D 2 
1 



-(iv k 2 +k%P 2 )t 



-0 



(< 



-(it) fc2+fcl-D2+2fciA;2Z)2)t 



fB 7") 

h> fc 2 + fc 2 ^ + 2k x k 2 D 2 ^ V 
Rescaling the integration variable h L = qi/Vi and expanding the resulting expression to 
0(1 /y/i) gives 

F( qi /Vi,q 2 /Vi) = 



v q 2 



q-iDi 



(« 



iv \/t 
q 2 D 2 2q x D 2 



-11- ^ 2 - gl 2 [ e -(i«0<?2V / t+g^2+2 9l < ? 2D2) _ ^ 

The limit t — » oo leads to the expression 

«^// c (°» c (°)(^- + 3) x 

m an ' 



(B.8) 



(B.9) 
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x e {qlW 2 )D2 e {qi+q2?D 2 ^-qxq 2 D 2 _ 2 + q ^J_^ 

The case of finite L can be treated by replacing F(q 1 /\/t, q 2 /\ft) by 



G{ qi /Vt,q 2 /Vi)= (B.10) 
1 (^-CmmVt+qlD2) _ 

^ e -(iv q2\/t+q 2 : D2+2q 1 q2D2) _ ^ 



iv q 2 + qiD 2 /Vt 
c 



qi 92 L 2 /t 



iv q 2 + q 2 D 2 / Vt + 2q 1 q 2 D 2 / V t 
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Figure Bl. Ensemble ( ) and effective (- -) dispersion coefficients of an exactly 

stratified medium for a point source; m = 10, D\ = 0.1, ctq = 1 




Figure B2. Ensemble ( ) and effective dispersion coefficients of an exactly 

stratified medium for a point source (- ) and line sources of various length (L = 1 
( ), L = 10( ),L= 100( )) ; t d = 10, D 1 = D 2 = O.f , tjg = f 
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Figure B3. Ensemble ( ) and effective (- -) dispersion coefficients for a stratified 

medium with a flow component transverse to the strata; td = 1, t v = 80, D\ = 0.1, 
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Figure B4. Normalised fluctuations of the dispersion coefficient in an exactly 

stratified medium for a point source ( ) and line sources with L = 10 (- -), 

L = 100 ( ), L = 1000 ( ); t d = 10, D 1 = D 2 = 0.1, a\ = 1 



